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Abstract 

We develop a rigorous analytical Hamiltonian formalism adapted to 
the study of the motion of two planets in co-orbital resonance. By con¬ 
structing a complex domain of holomorphy for the planetary Hamilto¬ 
nian, we estimate the size of the transformation that maps this Hamil¬ 
tonian to its first order averaged over one of the fast angles. After 
having derived an integrable approximation of the averaged problem, 
we bound the distance between this integrable approximation and the 
averaged Hamiltonian. This finally allows to prove rigorous theorems 
on the behavior of co-orbital motions over a finite but large timescale. 


1 Introduction 


Averaging methods are common techniques to study the dynamics of Hamil¬ 
tonian systems in celestial mechanics. The first and most famous example 
of averaged Hamiltonian system is perhaps the secular planetary problem, 
where the Hamiltonian of the planetary problem is averaged over the mean 
longitudes of the planets. The secular equations of the planetary motion 
appear in Lagrange’s work on stability of the solar system (Lagrange, 1778) 
while the secular Hamiltonian appears in Delaunay’s memory about the the¬ 


ory of the Moon (Delaunay, 1860). Poincare (1892) gave an expression of 


the secular Hamiltonian of the planetary three body problem while Laskar 


and Robutel (1995) present an analytical method allowing to compute the 
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expansion of the planetary Hamiltonian, and in particular, to get a concise 
expression of its secular part. 

The secular Hamiltonian of the planetary problem can be obtained, up 
to a finite order of the planetary masses, by averaging over the planetary 
mean longitudes (construction of a resonant normal form), the symplectic 
transformation that maps the non-averaged Hamiltonian to the secular one 
being close to the identity. Precise estimates on the size of these transfor¬ 
mations are required in order to prove the existence of invariant tori using 
the KAM theory. These kinds of rigorous estimates have been established 
especially by Arnold (1963), Fejoz (2004) and Chierchia and Pinzari (2011). 

When two planets are in co-orbital resonance, and more generally for 
two planets in mean-motion resonance, the transformation leading to the 
secular Hamiltonian is no more close to the identity, even at first order. 
Consequently, the secular motion does not provide a good representation of 
the real planetary motion. In these cases, it is still possible to use averaging 
methods, but for two planets the Hamiltonian is generally averaged over one 
fast angle, that is one of the planetary mean longitudes. 

Many authors work with the resonant averaged problem. Some of them 


(e-g. 

Erdi, 

1977; 

Namouni 

1999 

Morais 

2001, 

Robutel and Pousse 

2013), 

while others prefer a numerical averaging (e.g. 

Nesvorny et ah, 

2002 

Giup- 


and consequently the ’’distance” between the secular and the ’’complete” 
solution is rigorously estimated. 

In this paper, we estimate this distance and give an upper bound of 
the time for which this difference remains small enough. For this purpose, 
we derive, in section [2j a complex domain of holomorphy for the planetary 
Hamiltonian which allows to compute quantitatively the size of the trans¬ 
formations and of the perturbations involved in our construction. 

All the computations about perturbation methods are derived in section 
[3] and constitute the main novelty of this paper. 

The topology of the averaged Hamiltonian is studied in section [4] where 
we show the existence of an invariant manifold, associated to the quasi¬ 
circular motions, which carry an integrable dynamic. 

Section [5] is devoted to the construction of an integrable approximation 
of the averaged Hamiltonian and its degree of accuracy in the vicinity of the 
invariant manifold considered in section [4j We also give the general form of 
the solutions of this integrable system. 

Finally, in section [6j we can combine the quantitative estimates given in 
section [3] and the bounds on the remainders between the averaged Hamilto- 
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nian and its integrable approximation. This allows to prove rigorous theo¬ 
rems on the behaviour of co-orbital motions over large timescales. 

The last section concerns the proof of the technical propositions and 
lemma used in our reasonings. 


2 Hamiltonian setting of the problem 

2.1 Canonical heliocentric coordinates 


We consider two planets of respective masses mi and m 2 orbiting a central 
body (Sun, or star) of mass mo dominant with respect to the planetary 
masses. As only co-orbital planets are considered, no planet is permanently 
farther from the central body than the other, so the heliocentric coordinate 


system seems to be the most adapted to this situation. Following Laskar 


and Robutel (1995), the Hamiltonian of the three-body problem reads 


%(Rj,r,-) = UKCRj^j 

'H K (Rj,r j )= ^ 


+ 'Hp(R j ,r j ) 

with 

lh \ 

and 

,2 ft W r j\\J 



%p(Rj, r i) 


Ri ■ R2 

m 0 


mim 2 

y 71 - 

n - r 2 


(i) 


where r^ is the heliocentric position of the planet j, [3j = mom j (mo + 
m.j ) -1 and Jlj = Q(mo + fhj'), Q being the gravitational constant. The conju¬ 
gated variable of rj, denoted by Ry, is the barycentric linear momentum of 
the body of index j. In this expression, T-Lk corresponds to the unperturbed 
Keplerian motion of the two planets, more precisely the motion of a mass 
/ 3j around a fixe center of mass mo + rhj , while Tip models the gravitational 
perturbations. 

The Hamiltonian T~L is analytical on the whole phase space (M 8 since we 
consider the planar problem) except on the manifold which corresponds to 
collisions between two planets: 

2) = j(Ri, ri, R 2 , V 2 ) G M 8 such that ri 7 ^ r 2 j . 

If we introduce the small parameter e given by 


e = Max 


/ m\ m2 
\mo ’ m 0 


( 2 ) 
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one can verify that the Keplerian term of the planetary Hamiltonian is 
of order e and that the other one is of order e 2 as long as the mutual 
distance between the planets remain large enough. This feature justifies a 
perturbative approach. 


Remark 1 In the sequel, we will not give explicit estimates of the constants 
independent of the small parameters involved in the problem in order to avoid 
cumbersome and meaningless expressions. We will just denote uniformly by 
M a positive constant chosen sufficiently large for our purpose but indepen¬ 
dent of the relevant quantities, we will try to specify this at each occurrence 
of such constant but sometimes it will be omitted. 


2.2 The rescaled heliocentric coordinates 

According to |l]), the momenta Rj and the Keplerian part of the Hamil¬ 
tonian are of order e while the perturbation is quadratic in e. In order to 
get more homogeneous quantities, it is convenient to rescale the planetary 
masses, the Hamiltonian and the canonical variables. 

First, we introduce new planetary masses and new reduced masses by 
the relations: 


mi 771-2 a Pj 

7711 = —, rn 2 = — ; Pj = — 
£ £ £ 


--—-— and Hj = Q(mo+£mj) for j €{1, 2}. 

77io T £frij 


Moreover, for j £ {1,2}, we rescale the impulsions by R ; = efj while the 
positions r j remain unchanged. Hence, we have made a conformal symplectic 
transformation T which changes the symplectic form with 


drj A drj = e 1 ^ clUj A drj 
j j 


and the Hamiltonian linked to the considered system becomes 

U{v 3 ,y 3 ) = e~ l/ H o T(fj,rj). 


This leads to the expression: 

= 'H K {r j ,r j 

16(1,2} 


)+E'H P (Y j ,Y j ) 



with 

and 




r l ' r 2 n 771,1771 2 

- h n- 

171 o || r i - r 2 


( 3 ) 
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the canonical variables fj and r j and the Keplerian part Hk being now of 
order one. 

The three body Hamiltonian T~L is defined on 

D = {(?i, ri, r 2 , r 2 ) G M 8 such that ri / r 2 } 
since the positions remain unchanged. 


2.3 The Poincare variables. 


In order to define a canonical coordinate system related to the elliptic el¬ 
ements (a,j,ej, Xj,Wj) (respectively the semi-major axis, the eccentricity, 
the mean longitude and the longitude of the pericenter of the planet j ), 
we use complex Poincare’s rectangular variables (Xj, Aj, Xj, — ixj)j e <i 2 } £ 
(IxRxCx <C) 2 : 


A j = and Xj = \f~Aj\J 1 — \Jl — exp (ivjj), 


(4) 


this coordinate system has the advantage to be regular when the eccentric¬ 
ities and the inclinations tend to zero. 

Consequently, we have the product of the analytic symplectic transfor¬ 
mations around the circular orbits (i.e. for a given constant cq > 0): 




vV j j Xj 7 %xj ^ 


C 4 


O'= 1,2) 


(5) 


with 




(Xj, Aj, Xj. —ixj) where 


Xj G T 
A j G R* + 


and Xj G C with 



which yields the new planetary Hamiltonian: 


H(Xj,Aj,Xj, -ixj ) = H k ( Ai, A 2 ) + H P (Xj,Aj,Xj, -ixj). 
H is analytic on the domain 21 C (T x 1 x C x C) 2 defined by: 


^(Xj,Aj,Xj,-iXj) G <£j for j G {1,2} / 


$i(Ai, Ai,xi, -ix i) \ ^ ^ 
< h 2 (A 2 ,A 2 ,x 2 ,-ix 2 ) J 
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2.4 The 1:1 Resonance 

The Keplerian part of the Hamiltonian expressed in terms of Poincare’s 
variables reads: 

'H K {r j ,r j ) = H k ( Ai,A 2 ) = - ( 6 ) 

JS {1,2} ^ 3 

Consequently, the 1:1 mean motion resonance, also called co-orbital reso¬ 
nance, is reached when (Ai, A 2 ) = (A^, A°) such that the two mean motions 
are the same (let us denoted by ui this frequency), that is: 


dH K 

Mi 


(A?,A°) 


dH K 

dAo 


(A?,A§) 


:= uj > 0 . 


(7) 


or: 


(A?) 3 (A°)3 ' 


( 8 ) 


In order to work in a neighborhood of this resonance, we introduce a new 
coordinate system with an affine unimodular transformation completed in a 
symplectic way vH(Ci, C 2 , Zi, Z 2 , xi, x 2 ) = (Ai, A 2 , Ai, A 2 , xi, x 2 ) with: 


1 -1 

0 1 


Zi 

Z 2 


Ai - A? 
A 2 - AO 


(9) 


which yields slow-fast angles. 

More precisely, in a neighborhood of the co-orbital resonance, the angular 
variables evolve at different rates: £2 is a ’’fast” angle with a frequency of 
order 1, £1 undergoes ’’semi-fast” variations at a frequency of order yfe (see 
Section 5.2), while the variables Xj related to the eccentricities are associated 
to the slow degrees of freedom evolving on a time scale of order e (secular 
variations). With this set of variables, the planetary Hamiltonian becomes: 


H((j, Zj , Xj , — ixj ) = H o TfC/, Z 3 , Xj, —ixj) 

= H k (Zi, Z 2 ) + H P (Cj, Zj , Xj , -ixj). 


( 10 ) 


For an arbitrary A > 0, we consider the following domain centred at the 
1:1 keplerian resonance defined in the (£, Z, x, —ix) variables: 

4 0) = {((Ci, 0,0,0), (C 2 ,0,0,0)) € (T x M x C x C) 2 such that |Ci| > A} 
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where |.| denotes the usual distance over the quotient space T = M/27rZ. 
Hence, the angular separation A over /C^ yields a minimal distance between 
the planets: 


5 :=2asin 



where 


is the lowest radius of the orbits. 


a :=min 


(A?) 2 (A§) 2 \ 
Pirn ' film ) 


The condition |£i| > A for (i £ T can also be considered with the real 
variable (j g]A,27t — A[ since there exists an unique real representative in 
this interval for an angle with a modulus lowered by A. Hence has the 
structure of a cylinder in (]A, 2ir — A[xM xCxC)x(TxMxCx C). 

For p > 0 and a > 0 small enough (this will be specified in the se¬ 
quel), our initial domain /C^ can be extended in a complex neighbourhood 
CT C C 8 of the following type: 


4 c L={((^).o.^w. 2 } ec8 / ( 


(Re(Ci), 0,0,0) 
(Re(C 2 ), 0,0,0) 


e/C 


(o) 


and \Zj\ < p, |£j| < y/pa , \pj\ < y/pa, |Im(£j)| < a 


actually 



C 1C 


(€) 

A,p,£7 


nst 


= {(O’ z i> 0’ »/j)je{i )2 } G with 7 lj 



In this setting, we can define a complex domain of holomorphy for the 
planetary Hamiltonian where it will be possible to estimate the size of the 
transformations and the functions involved in our construction of a 1:1 res¬ 
onant normal form. More specifically, for A>0, p > 0, c > 0 and for p > 0 
we will consider the compact: 


/C. 


p ■ 


r (C) 

^ A,pp,pcr 


and the supremum norm ||.||oo on the space of holomorphic functions over 
the compact 1C V which will be denoted ||.|| p . 

With Hk of class on the compact {( 21 , 2 : 2 ) £ C 2 / max(|zi|, \z 2 \) < 
p}, there exists a constant M > 0 which satisfies: 


V (pi,p 2 ) £ N 2 such that \p\ \ + |p 2 | < 3, one has: \\ d zfCz^ H ^\\ <M l 11 ) 

1 ’ 2 n 


7 



where we denote ||.||p the supremum norm ||.||oo on the space of holomorphic 
functions over the compact {(21,22) £ C 2 / max(|zi|, |^21) < p}- 

By the real-analyticity of the transformation in Poincare resonant action- 
angle variables T = ($1 o 'P, <P 2 o 'I'), there exist po > 0 and 00 > 0 small 
enough such that T can be extended into a holomorphic function over the 
compact aQ with: 


f j o f : 


t-(C) 

(^ji Z j^ji r lj)je{l,2} 


<C 4 

(b> r i) 


for j = 1,2 


and the differential of *P is bounded by a constant C > 0 with respect to 
the supremum norm ||.||A,p 0 ,<n> on the space of holomorphic functions over 
the compact 

Without loss of generality, we can assume that C > 1. 


Theorem 1 There exist constants po > 0, cro > 0 and M > M independent 
of e and A (or equivalently of 6) such that if we assume: 

P< Po, cr <a 0 , 0 < p < a < 1 and a < sin < 5 (12) 

then the following bounds are valid: 

< Inf /Ci (\Hp(<;j,Z 3 ,Zj,rij)\) and \\Hp\h < M & - 

where we denote ||.||i the supremum norm ||.||oo on the space of holomorphic 
functions over the compact KL\ = /C^ CT . 

Remark: As it was specified, from now on we will denote uniformly by M 
a positive constant independent of e, p , a and A (or equivalently of 5) high 
enough for our purpose. We will try to specify this at each occurrence of 
such constant but sometimes it will be omitted 


3 Hamiltonian perturbation theory 

In this paper, we only consider the averaged Hamiltonian at first order in 
the planetary masses. More precisely, we prove quantitatively that there 



exists a canonical transformation C which maps the original Hamiltonian H 
in 


H ' (Cj, Z j ,x j ,-ix j ) = H K (Z 1 ,Z 2 }\Hp(C v Zj , xj , -ixj)+H,((., Z, } , x j , -ix j ). 
In this expression, the Keplerian part Hp reads 


Hk(Z i, Z 2 ) = - 


PM 




2(A? + Z 1 )2 ^l-Z^Z^' 


(13) 


while the averaged perturbation with respect to the second angle (which 
corresponds to a time averaging along the periodic orbits on the torus at 
the origin for the Kepler problem in our resonant action-angle variables) is 
given by: 

npiCi’Zj’Zj’-iVj) =%a / H p(.C v C . 2 + ^t,z j ,x j ,-ix j )dt 

' °r2i r (14) 

= bj i>C 2 , z v x v -ixj)d c 2 

v 0 

The remainder ,Zj,Xj,—iXj) is a much smaller general perturbation 

whose size will be estimated in the next section. 


3.1 Hamiltonian perturbation theory 

Now, we specify our construction of the averaging transformation which will 
be the time-one map C = of the Hamiltonian flow generated by some 
auxiliary function 

Using the Poisson bracket: 

L x(f) = ix, /} = dix-d 2 f + d 3 x-dif - d 2 x-dif ~ d 4 x-d 3 f, 

we have C = exp (L x ) and H = exp (L X )(H). In the new variables, the 
Hamiltonian can be written: 

H' = H k + H P + {x, H k } + {x, Hpj + H'-H- { x , H}. 

With a generating function x comparable with Hp, the terms of order 1 in 
this expansion (with respect to Hp) are : 

[Hp + {x, H k }\((., Z^Xj-iXj) = H P (CZ p x v -ix j ) 

+ VH K (Z 1 ,Z 2 ).^(( j ,Z j ,x j -ix j ) 
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(C) 

We will choose x as a solution, over /Ci = ’ Q , of the equation: 


dx _ _ _ _ 

=Hp(C 1 ,Z j ,x j ,-ix j ) - Zj, Xj, —ixj) 

actually, this equation is satisfied by : 


LO 


x(C j, Zj,Xj, -iXj ) = — J t. Hp(C 1 ,Z j , Xj, - iXj ) 


-#p(Ci ,C 2 + vt,Zj,Xj, —ixj) 


dt 


and 


H\Cj,Zj,Xj, —iXj) = 

Hk(? 1 ,Z 2 ) +Hp(C v Zj, Xj , -ixj ) + #*(£., Zj , Xj , -ixj ) 


with if* = (\7H k -7$) + {x, i/p} + H' — H — {\, if} for at : = 

9C 


3.2 Quantitative Hamiltonian perturbation theory 


Here, we make the construction described in the previous section with ac¬ 
curate estimates on the size of the Hamiltonian and on the size of the trans¬ 


formations which are involved. As in the section 2.4 for A > 0, p > 0 and 
a > 0, we consider the compact JC 2 and the supremum norm 11.|I2 on the 
space of holomorphic functions over /c2- With the expression of x and the 
fact that |[ffp|I2 < ||ifp|12? we obtain : 


2vr 

2 < —11 Hp 11 2 
a j 


(15) 


which allows to prove the following: 


Theorem 2 Let A > 0, p > 0 and a > 0 such that (A,2p,2cr) satisfy the 


condition (12), hence: 


p < f ; < ,<ir<mia (j’4c) for ,5;=sin (t) min (®? it) - (16) 
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In order to get a small enough transformation (or equivalently a small 
enough \), we assume moreover that 


£ (jJ 

bpa < 32-jtM 

Then there exists a canonical transformation C : IC 3 —> IC 2 such that 


(17) 


C is one—to—one and 1C 5 C CUCs) C 1C 7 (18) 

4 2 4 


and, still using the notations 

(0, Zj, Xj, —ixj)j£ { 1 , 2 } ~ ^ (Cjj Tp — > 
t/ie transformed Hamiltonian H' = H o C can be written as 


hi' (Cj, Zj , ) = -ffif {Z 1 , Z 2 /H p( ( v Zj, Xj , —ixj )+H*(C.., ^, -ixj) 

with || 77 p ||2 < 2 Mp and : 

d 

1177*11 3/2 < with 77 = 40^^ f- + (19) 

' 0 uj \a pad J 


The proof of this theorem is given in section 7.1 


The decreasing factor p in the averaged perturbation can be minimized 
for a fixed lower bound on the mutual distance <5 > 0 and a mass ratio s. We 


first relate the analyticity width p to S and e by choosing p = \j - such that 


the two terms in the factor p are of the same order, then 77 = 80 


Mn /e 
w ay/6 


For 6 < 16 ( 7 ( 70 , the maximal analyticity width a = 
factor: 


32C 


gives the minimal 


min(?y) = 2560 


CMn 


L) 


—7 for p = \ - and a = 


32C 


which imposes a lower bound 6 > (^(e 1 / 3 ) in order to get a decreased pertur¬ 
bation in the averaged system. Hence we recover the size of the Hill region 
inside which the averaged heliocentric Hamiltonian is not close to the initial 
one. 
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4 The averaged Hamiltonian and its topology 


4.1 Lagrange and Elder configurations in the averaged sys¬ 
tem 

Euler and Lagrange configurations are central configurations of the three 
body problenj^j The two Lagrange configurations correspond to the situ¬ 
ation where the three bodies occupy the vertices of an equilateral triangle 
(the equilibrium points L 4 and L 5 in the circular RTBP), while the three 
Euler’s ones are the aligned configurations (L\, L 2 and L 3 in the circular 
RTBP). In terms of heliocentric elliptic elements, these are represented by 
two homothetical ellipses. The motion of the planets on these fixed ellipses 
is Keplerian and the difference of their mean longitude £1 = Ai — A 2 is 
constant. More precisely, in the Lagrange’s case, we have: 

Ci = Ai — A 2 = — VJ 2 = ±7t/3, ei = e2, ai = 02 , (20) 

while Euler configurations lead to the relations 

Ci = Ai - A 2 = w\ - UJ 2 = 0, ei = e 2 , a± = a 2 + O^e 1 ^), (21) 

if the two planets are on the same side of the Sun [L\, L 2 ), and to 

Ci = Ai - A 2 = mi - zu 2 = 7r, e\ = e 2 , ai = a 2 + 0(e), (22) 

if the Sun is between the planets (L 3 ). In the both cases, these trajectories, 
in fixed reference frame, are periodic orbits whose period is the common 
mean motion of the planets. As a result, these central configurations corre¬ 
spond to fixed points of the averaged problem. More precisely, each type of 
configuration (Li to L 5 ) defines, in the averaged problem, a one-parameter 
family of equilibrium points. This implies that a given fixed point of one of 
these families possesses an eigenvector (tangent to the family) associated to 
a zero eigenvalue. This is the source of the degeneracy that will be discussed 
in section 15.31 

4.2 Some properties of the averaged Hamiltonian 

In this section we will study the main properties of the averaged Hamiltonian 
of order 1 : 

H(C V Zj,Xj, -ixj) = H'(Zj,Xj, -iXj ) - H*{Zj,^, -ixj) 

= h k (Z 1 , ) +Hp (C x , ?j, Xj , ~iXj ) 

1 See the webpages by A. Chenciner (2012) and by R. Moeckel (2014) at 
www.scholarpedia.org/article/Three_body_problem, and at 
www.scholarpedia.org/article/CentraLconfigurations^Euler. 
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and of its associated dynamics. 

The Hamiltonian Hp being an analytic function on the domain /C 3 , it 

2 

can be expanded in Taylor series in a neighborhood of (x 1 ,x 2 ) = (0,0) as: 
#p(Ci ,Zj,Xj,-iXj) = C'p.q {C 1 ,?i,?2)^ 1 vfx^ 1 xl 2 . (24) 

(p,p)sN 4 

In the previous summation, only the coefficients satisfying the relation 

Pl+P2=Pl+p2 (25) 


are different from zero. This propriety, known as D’Alembert rule, is equiv¬ 
alent to the fact that the quantity 


l?i | 2 + 


■f 21 


(26) 


is an integral of the averaged motion. The existence of this constant of 
motion makes possible the reduction of the averaged problem, leading to 
a Hamiltonian system which depends on two angles: the difference of the 
mean longitudes and the difference of the longitudes of the perihelion (see 
Giuppone et al.. 2010). This reduction, which decreases the number of 
degrees of freedom of the averaged problem from 3 to 2, introduces some 
technical issues (addition of a parameter, singularity when the eccentricities 
tend to zero). For this reason, we prefer not to reduce the problem. 


Besides the reduction, the relation (25) has many implications on the 
dynamics of the averaged system. 

One of these properties lies in the fact that the manifold 

Cq = {(C., Zj , aq, x 2 , x 1} x 2 ) G /C 3 , such that = x ■ = 0} (27) 


is an invariant manifold of the averaged Hamiltonian (23). Indeed, the 


relation (25) implies that the Taylor series (24) starts at degree two. As a 


consequence 


dxH P (Q , Zp 0,0) = dx Hp {C , Zp 0,0) = 0, 


(28) 


which proves the invariance of the manifold Cq by the flow of (23). 

Actually, with the previous definition, Cq is a complex manifold and we 
will consider the real manifold 


Cq = Cq n (] A, 27 t — A[xM xCxC)x(IxKxCxC) 


(29) 


which is also invariant by the flow of (23) since it is a real Hamiltonian. 
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4.3 The invariant manifold Co 
4.3.1 Hamiltonian dynamics on C o 


The dynamics on Co is given by the restriction of the averaged Hamiltonian 


(23) to this manifold, that is: 


H 0 (C V ?i, ? 2 ) — Hr(Z 1 , Z 2 ) Z 1 , Z 2 , 0, 0, 0, 0) (30) 


H'o =H 0 +H* 0 , 


(31) 


with 


TJ PlIH „ l COS C, 

no = —---- \-eymim2 

Zd i Za 2 


V a i a 2 a\ + a 2 — 2 aia ,2 cos ^ 


(32) 


and 


77 * Gmim 2 ( Pifa y/JnP 2 1 

fln = £ - --1 

•^0102 V m l m 2 G m 0 


cosC r 


(33) 


In the expression (32), a,j can be easily expressed in term of Z ,■ using 
the expression: 


®i = /h Vi -2 (A? + ? 1 ) 2 , a 2 = y 2 l p 2 2 (A° + Z 2 - Zy (34) 

deduced from Q. 

With the expressions [ij = rri 3 {l + 0(e)) (resp. Hj = Gmo + 0(e)) and 
the analyticity of the functions f(x, y) = x.y (resp. g(x, y) = y/x.y) over M 2 
(resp. M* xR*), we can write 


P1P2 V/^1^2 _ 
mim 2 Qmo 


< Me 


where M > 0 is independent of the small parameters in the problem. 


Moreover, we have a uniform upper bound on 


on 


cosC, 

A 1 A 2 


cos^ 


/ d\0,2 


(or equivalently 


over the compact /C 3 and gathering these estimates yield: 

2 


77; 


< Me 2 


(35) 
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4.3.2 Phase portrait on Co 

The phase portrait of the Hamiltonian Hq is represented on Fig. [l] in the 
plan (£ , u) where u is a dimensionless quantity related to the action Z x and 
to the semi-major axes by 

2w 1/3 /x“ 2/3 m\y/a\ — m2y/a2 x 

u= -— 1 - Z x = - - -v/^o + Oe) 36 

m i + m 2 2 

with /j ,o = QrriQ. This figure has been obtained for particular values of the 
masses and of the parameters given in the caption of the figure [lj but the 
qualitative structure of the phase space does not depend on these values. 

The phase portrait of i?o is invariant by the symmetry with respect to the 
Z 1 axis (Z 1 , £ ) i —> (Z 1 ,2tt — ( 1 ). When the two planetary masses are equal, 
the Hamiltonian Ho, and consequently its phase portrait, is also invariant 
by the symmetry (Z 1 ,C 1 ) > — > (-Z^CJ. 

The shaded areas indicate the outside of the co-orbital resonance. In the 
upper grey region where Z t > 0, the angle £ circulates clockwise while its 
circulation is anti-clockwise in the lower grey region (Z 1 < 0). The others 
domains correspond to the different kind of resonant motions. 

The two elliptic fixed points in the middle of two green areas, located at 


> r f /, r /,! , *7 £ mi — m 2 mim 2 2/3 - 1/3 , 2 \ 

C, G {vr/3, 571-/3} and Z 1 = - --- ^ u + 0[e ), (37) 

■ 1 b mi + m-2 mo 


are associated to the Lagrange equilateral configurations (see section 4.1). 
These points correspond to the relative equilibria where the three bodies 
occupy the vertices of an equilateral triangle rotation at the constant an¬ 
gular velocity u. Each of these points, named L 4 and L 5 in the RTBP, is 
surrounded by tadpole orbits (green regions) corresponding to periodic de¬ 
formations of the equilateral triangle. These regions, whose maxim al width 

for more 


Robutel and Pousse 


2013 


in the Z^^-direction is of order e 1 / 2 (see 
details) are bounded by the separatrix £3 that originates from the hyperbolic 
fixed point L 3 at 


Ci = 


£ mi - m 2 mim 2 2/3 -1/3 . w 

^ ?i = n -;- N u + 0 {£), 

2 m 1 + m 2 mo 


(38) 


for which the three bodies are aligned and the Sun is between the two 
planets and its separatrix. Outside this curve, in the blue domain, one 
find the horseshoe orbits that surround the three fixed points mentioned 
above. Along these orbits, the angle £ undergoes large variation such that 
£ [Cm, 27 r - Cm] with 0 < 0 (£ 1 ^) < Cm < 2 arcsin ((\/2 - l)/ 2 ) + 0(e). 


15 










The horseshoe region (in blue), whose vertical extend is of order e 1 / 3 , 
is enclosed by the separatrices associated to the two last Euler configura¬ 
tions located at £ = 0 ,Z 1 = ±0(e 1 / 3 ). These points, corresponding to 
the equilibria point L\ en L 2 in the RTBP, are associated with the Euler 
configurations for which the two planets are on the same side of the Sun. 

The last domain (in red), centred at the singularity of Ho, that is the 
collision between the two planets, is surrounded by the separatrix connecting 
the fixed point located at £ = 0 and Z 1 > 0 to itself. Inside this small 
region, the two planets seem to be subjected to a prograde satellite-like 
motion, the one revolving the other one clockwise. This last region (in red) 



Figure 1: Phase portrait of the Hamiltonian Hq in the coordinates 
The units and the parameter are chosen such that Z 2 = 0, Q = 1, u = 2 n, 
em\ = 10 3 , £ 777,2 = 3 x 1CF 4 . See the text for more details. 


and its neighbourhood that includes L\ and L 2 are located at a distance to 
the collision of order e 1 / 3 and consequently, is outside the validity domain of 
the resonant normal form (see Theorem[2]in the section 3.2) . Indeed, in this 
region the remainder H * is at least of the same order that the perturbation 
H P . 

Finally, let us remark that this figure is similar to the well known Hill’s 
diagram (or zero-velocity curves) of the non averaged planar circular RTBP 
(see Szebehely 1967) although the zero-velocity curves are not orbits of the 
system. It is also topologically equivalent to the phase space of the averaged 
planar circular RTBP when the eccentricity of the massless body is equal to 


zero (Nesvorny et ah, 2002 Morbidelli 2002) 
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5 An integrable approximation of the averaged Hamil¬ 
tonian 

5.1 Expansion around the resonance 

In order to get a more tractable expression of the averaged Hamiltonian Ho 
in the domaine /C 3 , it will be expanded in the neighborhood of Z ■ = 0 so 
that the remainder is always smaller than Hp. We will prove that the last 
condition is fulfilled when the expansion of the Keplerian part is truncated 
at the second order andH at zero order. 


5.1.1 Approximation of the Keplerian part 

Let us start with the Keplerian part of the Hamiltonian. If its constant part 
2 _1 (/3puf + /? 2 M 2 ^ 3 ) u2//3 i s omitted, Hp can be written as: 

Hk(Z 1 , Z 2 ) = wZ 2 + Q(Z l , Z 2 ) + Ri(Z 1 , Z 2 ) + R,2(Z 1 , Z 2 ) (39) 


( Ri(Z l ,Z 2 ) = Q(Z 1 ,Z 2 )-Q(Z 1 ,Z 2 ) 
with < 

1 R2^Z 1 , Z 2 ) = Hk{Z 1 , Z 2 ) — uiZ 2 — Q(Z±, Z 2 ) 
where the quadratic form Q reads: 

q(z 1 ,z 2 ) = -^ 4/3 (/3rvr 2/3 +/3 2 -v 2 ' 2/3 )zi 

+3u 4 / 3 P 2 1 f i 2 2/3 Z 1 Z 2 - ^ 3 P 2 l n~ 2/3 Zl 


(40) 


(41) 


and its approximation Q(Z 1 , Z 2 ): 

Q(ZpZ 2 ) = -^ 4 / 3 / x 0 - 2/3 (mr 1 +m 2 1 )Z 2 1 

+3u 4/3 HQ 2/3 m 2 1 Z 1 Z 2 - ‘^ A/3 iiQ 2/2 ‘m 2 1 Zl, 


(42) 


We first look at the terms R,\ , as: 

/ 3j = rrij + O(e) and nj = Qmo + 0(e) = no + O(e), 
the quantity R\ = Q — Q satisfies the relation: 

R 1 = 0(e\\{Z 1 ,Z 2 )\\ 2 ) (43) 
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and more specifically, there exists a large enough constant M > 0 indepen¬ 
dent of s, p, a and 6 such that 

||fli||3 < Mep 2 . 

2 


As regards the estimate of i? 2 , the application of the Taylor formula on 
the function g(t ) = Hx{tZ 1 , tZ 2 ) for (Z 1 , Z 2 ) £ AA leads to: 


Z 2 ) — Hk(Z 1 , Z 2 ) — ujZ 2 — 


Q(Z 1 ,Z 2 )= f 

■Jo 


1 (1 - *) 2 


g {3 \t)dt 


which finally leads to: 


Using the inequality (11), we have for all t 6 [0,1]: |<p)(t)| < M\\(Z 1 , Z 2 )|| 3 , 

R 2 (Z 1 ,Z 2 ) = 0 3 (Z 1 ,Z 2 ), (44) 

9 M 


and more specifically to: II-R 2 1 |a < - P 3 ■ 

2 16 

5.1.2 Approximation of the perturbation 

We consider the function G{C , x-,x) on the compact JC 3 which is equal 
to the modified function Hp(£ , 0,0, Xp Xj) where the resonant actions A° 

equal to /x^ 3 /3,-u; -1 / 3 are replaced by g?Q 3 f3jUJ~ 1 / 3 . Equivalently, using the 
notation: 


[go 3 / Vj 


3 /^ 


AZ“ = .6 hT W 12.0 _ W 12.0. = 0( e) for j e {1, 2} 


CJ 


yUj V W 

and using the transformation ([ 9 ]) we obtain: 

G(C 1S =f?p(C 1 , AZ°, AZ° + AZ 2 ,Xj,Xj). 


(45) 


Then the averaged perturbation 2fp can be split in the sum of three term 
as follows: 

JJJ P (Cp ? li j i "P j ) [ J (C-^ 1 P j 1 ) T ^3 (C-j^ 1 Pj 1 ) T -^4 (C-^ ? ?j 1 J' j 1 j ) 

(46) 

with 


^3 (C 1? ji 5 ) -^p (Ci? o ? 0 ^ (C 1 ? *? j ? j■ ) 

A 4 (C,: Z^X^Xj) = #p(Cj, ^ 1 , Z 2 , Xj , Xj) - Hp (^ , 0, 0, X j , Xj) 


(47) 
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In order to estimate theses remainders, we need to consider the smaller 

(k) 

compact /Q' C /C 3 where the neglected terms are smaller than the averaged 

2 2 

perturbation I i=> 113 ^ < Il-Hplk- 

2 _ 2 

With similar reasonings as in the previous section, we use the mean value 
theorem to evaluate the remainder in the truncation at order 0 of H p we 

(k) 

consider ICp C JC P which is defined for p > 0 and k > 0 by: 

i )2 } G fcp such that max(|Zi|, \Z 2 \) < Kpj . 

The supremum norm ||.||oo on the space of holomorphic functions over the 
compact will be denoted by ||.||p K \ For p = 3/2, we obtain on this 

smaller compact K.\ ; C /G for a small enough k > 0: 

2 2 


l|i? 3 |lS K) = I 0,0, Xj,Xj) -Hp(C v AZ^, t\Zl + AZ%,x j , Xj )1£ 


7 ^ ui(«) 


which implies that 

\m\i K) < WdzHplkmz^AZ^ + AZ^W < 2 |^p|| 2 ||(AZ 1 0 ,AZ 1 0 + AZ 2 °)|| 

2 2 p 

We have an upper bound IKAZj 5 , AZf + AZ 2 )|| < Ce for some constant 
C > 0 and our estimate on Hp yields: 


||ffp|| 2 <2M§=> \\R 3 \\i K) <M— (48) 

0 2 pO 

for a large enough constant M > 0. 

In the same way, the mean value theorem yields: 

m\\i K) < 11^11311(^,^)11 < 2 |^p|| 2 ||(Z 1 ,Z 2 )|| 

2 2 p 


11 -^ 4 (C/ ? ^\ 1 4^2 1 5 4 / 


3 ^ < AM^—np = AaPk. 

2 pO 0 


(49) 
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5.1.3 The final approximation of H 


Gathering the approximation given in sections 5.1.1 and 5.1.2, and omitting 
the constant term 2~ x + /? 2 Aw 2//3 , the averaged Hamiltonian if 

takes the following form: 


= u;Z 2 + Q(Z 1 ,Z 2 ) + G{C v x j ,x j ) + R{C Zj, x j ,Xj), (50) 


where the quadratic part Q is given by the expression (42), the perturbation 


G is defined in the top of Section |5.1.2t and the remainder R is defined by 
the sum R = R\ + f ?2 + R 3 + R 4 where the Rj are given in (40) and (47). 
If, as in the end of the section 3.2, we relate the analyticity width p to 6 


and e by p 2 5 = e, we get the following upper bound with M large enough: 



since p < 1 and 5 < 1 =>■ e < p 2 . 

Especially, we have Il-Rll^ < 3Mp 3 + AM up 2 < mp 2 < ||ffp||3 < 

2 2 

|(ffp|| 3 K ) f° r K small enough and e small enough (since p < e). 

2 


5.2 The dynamics on C 0 and its implication for the initial 
problem 


Using the approximation (50) of the averaged Hamiltonian where the re¬ 
mainder R((jZj,Xj,Xj) is neglected, its restriction to the invariant manifold 
Co reads: 

R 0 = ojZ 2 + Q(Z 1 , Z 2 ) + (52) 


m 0 


with 


F(C 1 ) = G(C 1 ,0,0) = cosC 1 



(53) 


In order to uncouple the fast and semi-fast degrees of freedom, we define 
the symplectic linear map (<+ £ , Z 1 , Z 2 ) = £(<pi, ip 2 , h, h) on the cylinder 
] A, 27 t — A[xT xRxl by: 



1 

mi 

mi+7712 




(54) 


20 
















and completed by the identity in the Xj and x 3 variables. As a consequence, 
we have: 




_ ,4/3. -2/3 

~~2 


+ uh - -U) ' 11 o 


, 4 /3 2/3 


1 1 

-1- 

mi m 2 


j2 . 2/3 2/3 m l m 2 z7 , x 

/l +^o / - / —^i) ( 55 ) 


mi + m 2 


The dynamics of the fast variables (</> 2 , H) is now governed by Hq ; while 

the dynamics of the semi-fast variables (</q,Ii) is given by Hq 1 *. Of course, 

( 2 ) 

these two Hamiltonians are integrable. But the dynamics of Hq is trivial 
while that of is less. 

In the domain K,^\ the phase portrait of the two approximations 


(already explored by several authors, i.e. Morais, 2001, and references therein) 
and Ho of the averaged Hamiltonian are topologically equivalent. They both 
have two elliptic fixed points corresponding to L 4 and L 5 and an unstable 
equilibrium associated to the Euler configuration L3. The stable equilib¬ 
ria are located at = (±vr/3,0), and their eigenvalues, which are the 

same for both points, are equal to mi +™ 2 uj+O( e 3 / 2 ). The unstable 

point is located at (</?i,ii) = (715 0), its eigenvalues are =fcy ^ 1711 w + 
0(e 3//2 ). As we have seen in the section 4.2, in the domain enclosed by 


the separatrices emanating from this fixed point, one find tadpole orbits 
surrounding L 4 or L 5 . Outside these invariant manifolds are the horseshoe 
orbits that encompass the three fixed points mentioned above. 

At this points, we know at least qualitatively what are the orbits on 
the invariant manifold Co- To go further, we would like to have the tem¬ 
poral parametrization of the corresponding trajectories. However, even if 
the Hamiltonian T~L^ is integrable, its trajectories cannot be given explic¬ 
itly. Consequently, in the sequel, we will assume that these trajectories that 
satisfy the canonical differential equations: 


h 

<pi 


= e/ig^w 2 / 3 — 1 —- (l — (2 — 2cos^i) 3 / 2 ) sin</>i 
m 0 V / 


= 3 mi + m 2 , . 4 / 3 j) -2/3 


mi m 2 


0 ' h 


(56) 


are perfectly known, once given its initial conditions (</>i(0), /i(0)). Ac¬ 
tually, these solutions are all periodic since the level curves of Hg 1 "* are closed 
and without singularities. 
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5.3 The normal stability of the manifold Co 


Now, we study the linearized dynamic around the invariant manifold Co- 
Consider an arbitrary trajectory on Co- It is shown in Robutel and Pousse 


(2013) that the variational equation in the ( Xj,Xj ) direction around this 
solution reads: 

X = M(t)X 


where 


X = 


and M (t) = ieuj 


m\m2 

mo 


Mcy)) 

m\ Jm\m 2 
\ y/mini2 m 2 


with 


4D( C x )5 

DiQ = y/' 2-2cosC l5 


^5 cos 2(' | —13 + 8 cos — cos 

(e“ 3 ^i + 16e“ 2 ^i - 26e"^i + 9e*0 


(57) 


(58) 


(59) 


£ (t) being a solution of the canonical equation associated to the Hamilto¬ 
nian (52) 

According to the Floquet theorem (see Meyer and Hall, 1992), if the 
frequency of the considered periodic solution is v, the solutions of the vari¬ 
ational equation take the form 


Y(t ) = P{vt) exp(Ht), 


(60) 


where U is a constant matrix and P(tp) is a matrix whose coefficients are 
2vr-periodic functions of ip. As, if Y is a fundamental matrix solution to the 
variational equation along a 27r/z'-periodic solution, one has the relation 

Y(t + 27rzz -1 ) = Y(t) exp (2^ 1 I7) . (61) 


It turns out that the stability of the solutions of the variational equation 


(57) depends on the eigenvalues of the matrix U. As stated in section 4.2 


the quantity x 1 x 1 + x 2 x 2 is an integral of the variational equation (57). This 
implies that the solutions of (57) are bounded, and as a consequence, U is 
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diagonalisable and the real parts of its eigenvalues are equal to zero. Al¬ 
though we cannot exclude that one of the eigenvalue vanishes, the invariant 
manifold Co is normally stable. 

There exists, on Co, at least three trajectories for which one eigenvalue 


of the variational equation (57) vanishes. These are the stationary solutions 


of the differential system (56), that is the Euler configuration L 3 and the 
Lagrange ones L 4 and L 5 . 

For the collinear configuration associated to L 3 , which corresponds to 
= (-7T, 0), we have 


7 


^ 0 ) = B(n) = g. (62) 

As a consequence, the equilibrium has two eigendirections collinear to 

(63) 


V l = ( V™2 


y/m{ 




associated respectively to the eigenvalue 
1 . 7 mi + m 2 

v n = a- 

8 mo 


and vZ = 0 . 


(64) 


The reason why one of the eigenvalues vanishes has been given in section 
4.1 Indeed it is easy to verify that the eigenvector corresponds to the 


elliptic Euler’s configurations where the two planet are in the two sides of 
the Sun. 

The other eigendirection corresponds to a non trivial family of periodic 
orbits. According to the expressions ( f63| ), the configurations corresponding 
to Vj are two ellipses in conjonction (zu\ = ro 2 ) with equal semi-major axis 
and whose eccentricities satisfy the relation m-iei = m 2 e 2 . 

Contrary to the previous Euler’s configurations, the ellipses are not fixed, 
but precess at the same rate defined by the frequency —iv\ = 0(e). This 
eigendirection gives rise to a one-parameter family unstable periodic orbits of 
the averaged problem (periodic in rotating frame in the non-averaged prob- 


lem) described by 

Hadjidemetriou et al. 

(2009) and related to the Poincare 

solutions of second sort (see 

Robutel and Pousse 

2013 

for more details). 


Similar phenomena occur in the neighborhood of the two equilateral fixed 
points L 4 and L 5 at ipi = ±7t/ 3 and I\ = 0. Without entering into details 


(see Robutel and Pousse, 2013), let us just mention what we get for L 4 (the 
results are similar for the other equilateral equilibrium). The coefficients of 


the matrix (58) satisfy 


A(—) = - — , 

V 8 ’ 


27, 


^ 3 )=1 6 ( 1 -^)- 


(65) 
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As a consequence, its eigenvectors are 


v i = ( 

V — \A«i 


and V?,,= 

n ' 3 \ y/nV2 


ii r/3 


and the corresponding eigenvalues are 

27 mi + m 2 


t/3 


= — l£UJ 


8 


and 


m 0 


t/3 


= 0 . 


( 66 ) 


(67) 


It is easy to verify that the configurations associated to CC are the elliptic 
equilateral configurations that are fixed points of the averaged Hamiltonian. 
This is the reason why = 0. As in the case of the Euler configuration, 
the eigenvector is tangent to a one-parameter family of periodic orbits, 
called anti-Lagrange by Giuppone et al. (2010). For small eccentricities, 


the elliptic elements of the corresponding orbits that precess simultaneously 
at the frequency iv ] n / :i satisfy the relation miei = m 2 e 2 and w\ — V 02 = 
Ai - A 2 + 7T. 


6 Consequences for the co-orbital motion 

So far, the study of the manifold Co and its dynamics was carried out in the 
context of the averaged problem. To end this section, we will study what 
happens to the dynamics of Co when we go back to the initial variables. The 
following theorem gives a partial approximation, on a finite but large time, 
of the dynamics on the manifold Co in the initial variables —ixj ). 

Actually, we explain at the end of this section the way to obtain a complete 
result of approximation with our methods. 

We consider a solution in the initial variables (\j(t),Aj(t),Xj(t), —ixj(t )) 
starting with at t = 0 in the image of the manifold Co by the transformation 
(\j,Aj,Xj, —ixj) = T o C o Ij,x j: —iXj) considered in section 2, 3 and 
5, hence: 


(Aj(0), Aj(0),Xj(0), -ixj( 0 )) = $oCo%(0),/ i (0),0,0). 
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Using these linear and averaging transformations, we can formally write: 

Ai (t) = u(I 2 {0))t H-^— ipi(t) + p 2 (0) + pi(t) 

m i + 7712 

\ 2 (t) = u(I 2 (0))t - ^ - P! (t) + p 2 ( 0) + p 2 (t) 

777 1 + 777-2 

Ai(t) = A? + h(t) + m2 I 2 { 0) + n(t) (68) 

mi + m 2 

A 2 (t) = A° - /!(t) + mi / 2 (0) + r 2 (t) 

7?7l + 777-2 

*j(<) = w(t) 


where the function (p\{t), I\(t)) in these expressions is the solution of the 
differential system (56) with the initial conditions given as follow: 


¥>i(0) = Ai(0)-A 2 (0) 
777 2 


h( 0 ) = 

¥>2(0) = 


7771 + m 2 
7771 


(Ai(0) - Aj) - 


mi 


Ai(0) + 


m 2 


mi + m 2 

A 2 (0) 


(A 2 (0) - A") 


mi + 7772 mi + 777-2 

/ 2 (o) = a 1 (o)-a; + a 2 (o)-a° 

w(/ 2 ) = d l2 n ( 0 2) (p 2l I 2 )=u- 3w 4/3 /7o 2/3 


h 


7771 + 7?7 2 


(69) 


and the remainders pi(t), p 2 (t), T\ (t). r 2 (t), Tj+ 2 (t ) should be small over 
large times. 

We give a partial theorem in this direction. 


We first recall bounds on the remainders in the previous computations 
with the choices 5p 2 = e and 326*0- = 6 < 16600 , there exists a large enough 
constant M independent of e, p, a and 5: 


and we denote 


\H*\\3/2<M X / T5 


m 


= eO(e) = 0(e 2 ) < Me 2 


\\R\\i K) <Mp (3/7 2 + 4kp) 


(70) 

(71) 

(72) 


H* +H* 0 + R 


(«) 

3 

2 


<L:=M 



+ £ 2 + p (3p 2 + 4k/5) 


(73) 
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Theorem 3 We denote: 


mi + m2 
2m\ + m 2 


and T p the first time of escape out of 


for the solution (</?,,■ (0), lj(0), 0,0) along the flow governed by the av¬ 

eraged Hamiltonian H'. 

With the previous assumptions and notations, we have the following 
bounds on the remainders: 

\r j+ 2 {t)\ < p for |t| < min (74) 


and 


Ai (t) + A 2 (t) = Ai(0) + A 2 (0) + r(t) 


with |r(t)| < p for |i| < min ( T p 


np 


(75) 


The proof of this theorem is given in the section 7.3 and exactly the 
same reasonings would give a complete approximation of the solutions on 
Co except that we need moreover the action-angles variables linked to the 
integrable Hamiltonian Hq 1 which can be built by classical techniques (cf 
Arnold) 

The first step is to bound the difference between the flows linked to two 
nearby Hamiltonian in the normalized variables (jpj, —iXj) £ fcp'' 1 and 

then gives a sharp timescale such that these two flows remain close in the 
initial variables (Xj, Aj 


7 Proof of the theorems 


7.1 Proof of theorem [l] 

Using the notations of section [2~i] and the mean value theorem, we obtain a 
lower bound for the minimal distance (with C 2 equipped with the Hermitian 
norm) between two orbits in the complex domain IC A po aQ . 

Actually, denoting 


(?i 


(C) (C) -(C) (C) 


■ 1 


) = T tej,Zj,£j,rij)j e{1 ' 2} with T = ($1 o 1> 2 O 'k), 


26 




we have: 


= *((Re( ? i),0,0,0);(Refe), 0,0,0)), 


and the mean value theorem together with our bound C on the differential 
of T allows to write for j £ {1,2}: 


I \(Sj,Zj,€j,r}j) ~ (Refe),0,0,0)|| < p + 2y/pa + a 
=$■ r j C) -rj R) < C(p+ 2^/pa + a) 

then the triangular inequality yields: 


(C) (C) 


r, — r 


> 

> 


r (C) r (R) , r (R) P (K) , r (R) _(C) 

r l r l ■+■ 1 1 I 2 “T r 2 1 2 


r - r 0 


r - r.; 


_(C) «) 


_(C) _(R) 

r i — r 1 


r 2 “ r 2 


— 2C(p + 2y/ per + c ) 


(76) 


If we assume that 


6 


p < aand ( 7< —= -sin(-), 


A 


(77) 


we obtain with the lower bound on the mutual distance in the real domain: 


„(C) „(€) 


<5 


> 6 — 8Ca > - = a sin — 


A 


and the expression of the planetary Hamiltonian gives the upper bounds on 
the size of Hp in the complex domain. 

In the same way, under the assumption p < a < 1, we obtain with the 
upper bound on the mutual distance in the real domain: 


„(C) JC) 


- r.. 


< M + 8 Ca < M + 8C 


(78) 


and the expression of the planetary Hamiltonian gives the lower bounds on 

the size of Hp in the complex domain provided that M is large enough to 

satisfy — < M + 8 C. 

M 
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7.2 Proof of theorem [2] 

Let A > 0, p > 0 and a > 0 such that (A, 2p, 2a) satisfy (12), hence: 


p < a < 


32 C 16C 


A 


sin — with 5 := 2a sin — 


A 


(79) 


then there exists a constant M > 0 independant of e, A, p and a such that 


\\Hk\\ 2 P < M and \\H P \\ 2 = ||/7p||a,2 P ,2 ct < M- 


(80) 


for the supremum norm ||.|| 2 P (resp. ||.||2 = ||.||a,2 P ,2ct) on the space of 
holomorphic functions over the compact 

{(Zi, Z 2 ) G C 2 such that max(|Zi|, \Z 2 \) < 2p} 

(resp. over the compact 2a = K 2 ). 

In view of our bounds (15) and (80), we obtain: 


2ir 2 tt£ 

\\ x \\ 2< —(\\H p \\ 2 )<M— . 

0J 0J0 


(81) 


Moreover, with our rescaling of the masses, the derivatives of the Ke- 
plerian part Hk remains of order one and we can also assume that the 
constant M > 0 independant of e, A, p and a gives also an upper bound on 
the Hessian of Hk : 


\/(Zi,Z 2 ) € C 2 with max(|Zi|, \Z 2 \) < 2p 


d 2 H K 

dZidZj 


(Zi,Z 2 ) 


< M. (82) 


We must first estimate the size of the partial derivatives and the Poisson 
bracket of analytical functions on /C 2 by classical applications of Cauchy 
inequalities which can be found in Poschel (1993) and Giorgilli (2003). 

Let / be analytical on /C 2 (continuous on the boundary), we can write : 


d£ 

d( 


(C> Z, x, —ix) 


Su P||(ei,e 2 )|| = l 


d 

— /(C + te, Z, x, —ix) 

dt \t=* 


for (C ,Z,x,-ix) = (Cj, Z 3 px r/ G IC 2 . 
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One then applies Cauchy formula to the function 1 H > f(C + te, Z, x, — ix ) 
of the complex variable t, holomorphic and continuous on the boundary for 

/(P\ 

\t\ < a/4 when (£, Z, x, —ix) E JCy 4 , and obtains : 


0J_ 

9C 


< 


2 • 


7/4 


<7 


The same reasoning yields the equivalent inequalities for the other partial 
derivatives : 


Of 


dZ 


7/4 


4 

< - 
P 


Of 

and 

Of 

dx 

7/4 

dx 


< 


7/4 


sfpO 


To estimate the size of the Poisson brackets for two analytical functions 
on JC 2 , we write in a similar way that 


{f,g}((,Z,x,-ix) = 


<1 

dt\t=o 


. df df . df df 


and the function of the complex variable t in the right hand side is defined 
for \t\ < (we apply the Cauchy formula to /), hence we can write: 


||{/, 2 } 117/4 < 16 


II/II2IHI2 

pa 


in the same way we obtain 


I|{/j 1 3/2 < 4 



2 


Since the averaging transformation (p is the time-one map 4*/ of the 
Hamiltonian flow generated by some auxiliary hamiltonian for any func¬ 
tion K defined on the phase space: 


j t (K o <&*) = {K, x} o ^ => (. - 0 = £ ^ o 


Cj Cj 


3/2 


< 


d'X 


dZi 


7/4 


for j e {1,2}. Actually, we prove that starting inside /C 3 / 2 along the flow 
yields a time of escape out of /C 7/4 which is bigger than 1 and we can 
use our estimates on x- 
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Hence , we have: 


C, - Ci 


4 8Mtt e a 

< -llxlb <-i' < 7 

7/4 p UJ pd 4 


with our threshold ©• 

In the same way, we obtain: 


\\z- - zA\ , < 

dx 

W-3 -7 1 13/2 — 

d(j 


7/4 


1 ^-^’ 11 3/2 ^ 


l?i-^ll 3 /2 ^ 


d\ 


dxj 


dx 


< 


dXn 


7/4 \[PV 

4 


4 8Mn s p 

< |Ixl12 <- - < ~ 

G UJ GO 4 


8Mir £ J~pG 
2 <-— < 


< 

7/4 yfpv 


UJ y/pcfd 4 

8Mtt £ yfpd 

llxlh <-< — 

UJ pGO 4 


which yields JCs C C (1C 3 ) C JCt . 


As H'iCj, Z r /x 3 , -iXj) = H K (Z l , Z 2 ) 

+Hp(C 1 ,Z j ,x j , -ixj) + H, ((., Zj , Xj , -iXj ), 


where 

H* = (VH K -Tf) + { X , Hp} + H'-H- {*, H} for := ( ° ) , 


we may use Taylor’s formula at order two to write: 


H'-H-{ X ,H} = Ho$*-Ho$X- 



(1 -t){x,{x,H}}°®tdt- 


Since 4>^(/C 3 / 2 ) C /C 7 / 4 , using the upper bound M on the Hessian of Hp 
and Cauchy inequalities, one finds the following estimate : 


moreover : 


(X/Hk-^). 9 ^ 

d( 

I |{X) Hp} 11 7 /4 < 16 


7 | |xl I 2 77 tM 2 p£ 

< M-p——— < -— 

4 cr/4 uj gS 

7/4 ' 

||xl| 2 ||gp|| 2 < 16 vrM 2 e 2 

pG ~ uj pg 5 2 


(83) 


(84) 
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To estimate the third term in H *, we insert again the definition of x 
given in section (3.1) into the Poisson bracket to get : 

{x, H}=Hp-H P + ( V H K - + {x, H P } 

d( 

hence 

\\{X,{X,H}}\\ 3 / 2 <\\{ X ,Hp-Hp }}\\ 3/2 


+ 


x, (VHk-12).^ 


+ \\{x,{x,Hp }}\\ 3/2 


3/2 


Using again Cauchy inequalities, we can estimate the previous terms and 
the sum of these inequalities yields the given value of r), indeed: 

nr Tt Tin „ JlxlUWHp —Hp \\2 16itM 2 e 2 

\\{X,Hp -tfp}}|| 3 /2< 4 - < 


X,(VH K -^) d -± 


pa 

I Ixl | 2 


co pad 


2 ' 


< 


d( 


7/4 


3/2 


pa 


< 7 M 


47 tM pe 
co ad 


ll{x,{x^p}}||3/2<8 " /vl|2|l{X,gp}ll3/4 <M 


( 167 tM 


pa 


y oo pad I d 


moreover the threshold 


£ CO 

— < - P 

da 327 tM 


< 


CO 


dpa 32ttM 


gives: 


x, (vh k - rt). 


dx 

d( 


< 


7 7 tM 2 p £ 
2 co a d 


< 


00 


I |{X) {Xi Hp \} 113/2 < 


3/2 

87 tM 2 £ £ 


dpa 32irM ' co pad d 

Finally, the sum of all the previous estimates yields: 

11-^*113/2 < 


CO 


Mn f 21 p e \ e 
~2 a +4 °^d) M d 


and ensures the claimed value for the upper bound p = 40- 


(85) 

( 86 ) 


. Mn ( P £ \ 
= 4n- - a - 


00 \a ^ padJ 
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7.3 Proof of the theorem [3| 

The first part of the proof of theorem [3] comes from the use of Cauchy 
inequalities to bound the size of the Hamiltonian vector field linked to (77* + 
Hq + 7?) o C over the compact JCp K \ 

The second part comes from the size of ToCo£ composed of the averaging 
transformation C and two linear transformations T and C and the choice of 
a time T which give terms of the same order in the upper bound on the 
error in the approximate solution. 

More specifically, by straightforward computations the norm of the trans¬ 
formation C admits the upper bound p^ 1 = ^ ience 

C C JC[ k) . 

With our notations, the function 

(<7h(7), h(t), 7 2 (t), 0,0) = (</?i(t), ^2(0)+3/ 2 ^ 2) (/ 2 (0))t, h(t), 72(0), 0, 0) 

(yi(*), ¥>2(0) +(u>- 3u; 4 / 3 Mo 2/3 ^ ) t, h(t), 7 2 (0), 0,0) (87) 

v mi+ra 2 / 

( 2 ) 

is a solution of the Hamiltonian system linked to 77 q ' on Co, hence it is a 
solution of the averaged system governed by 77. 

Then, the complete averaged Hamiltonian H' satisfies 77' —77 = 77* + 

77q + R and if we denote h(t), 72(7), x(t), x(t)) the solution of 

the system linked to H' starting at (<£>i(0), <^ 2 (0), 7i(0), 72(0), 0,0), we can 
write _ 

d ^ = -d^+Hl + R)ojC. 

Then, with our assumption on the time of escape T p , we have: 
£(^i(7),^2(7),7i(t),7 2 (7),x(t),x(7)) eK,[ K ^ for |7| <T V 


and Cauchy inequality yields: 


dl 2 

2 


< — 

R* +T/q + R 

dt 

up 


(«) 2 L 

3 < —• 
I K P 


(88) 


— — 2 L 

h{t) = 7 2 (0) + r 2 {t) with |r 2 (i)| < — \t\ for |7| < T p . 

up 


We have the same bound in the averaged variable Z 2 = 7 2 . 
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Finally, in the initial variables Z 2 . we have: 


\?2 


^2113/2 - 


dx 

d(2 


< 


7/4 


8Mir e 
uj a 5 



et ||Z 2 (t) - Z 2 (0)|| < ||Z 2 (t) - Z 2 (t)\\ + \\Z 2 (t) - Z 2 (0)|| + ||Z 2 (0) - Z 2 (0)|| 
=> 11^2 (t) ~ Z 2 (0)|| < 2 ||Z 2 - Z 2 || 3/2 + I 2 {t) - hit)) 


=► 11^2 it) - z 2 (0)|| < ^ + ^|t| < p/or |t| < min (89) 

which gives the formula for A in the theorem. 

The formula for (xj, x-j) is proved in the same way. 
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